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Saturation  and  the  “Universal”  Spectrum  for  Vertical  Profiles 
of  Horizontal  Scalar  Winds  in  the  Atmosphere 

E  M.  Dewan  and  R.  E  Good 

Atmosfherie  Oyties  Branch,  Air  Force  Geophysics  Laboratory,  Hansom  Air  Font  Bat*.  Manachuuru 

—f-S  A  theory  it  presented  which  espial  ni  iho  universal  aalurt  of  ruw-dimemionil  vortical  wave  number,  it, 
power  spectral  densities  (PSDs)  of  horizontal  winds  as  measured  in  the  atmosphere  and  predicted  by 
VanZandl.  The  theory  is  that  the  PSD  amphiude  m  any  liven  wave  number  (greater  than  a  certain 
minimum.  *.)  is  determined  by  its  saturation  value  due  either  to  shear  instability  (ia,  critical  Richardson 
Number)  or,  more  likely,  to  convective  instability.  This  explains  why  the  PSD  amplitudes  observed  do 
not  grow  exponentially  with  increasing  altitude. s^his  saturation  theory  assumption  plus  other  consider¬ 
ations  leads  to  a  PSD  of  the  form  Nr/k“,  where  als-iqthe  range  of  about  25  to  3  and  N  is  the  Brunt 
frequency.  A  simplified  model  involving  superimposed  nafteaufepads  of  gravity  waves  as  well  as  a  model 
based  merely  oo  dimensional  arguments  both  lead  to  a  -  3?*Tbe  full  model  not  only  captains  the 
observed  spectral  slopes  but  also  predicts  the  PSD  amplitude  in  the  troposphere  to  be  3 5  times  smaller 
than  in  the  stratosphere.  )The  derivation  of  the  model  is  baaed  on  the  saturation  condition  that 
J  k’PSDlk)  4k  m  w*.Sf he  model  may  also  apply  to  the  ocean  and  explain  the  Garrett-Munk  vertical 


wave  number  spectrum. 


-  (jro //<") 


1.  Introduction 

In  previous  communications  {Dewan  ti  al ,  1984]  it  was 
shown  that  the  vertical  wave  number  power  spectral  densities 
(PSDs)  of  horizontal  winds  in  the  stratosphere  were  nearly 
identical  over  the  wavelength  range  of  40  m  to  1  km.  This 
occurred  in  spite  of  the  fact  that  the  data  were  obtained  at 
different  times  of  day  (dawn  versus  dusk),  during  different  sea¬ 
sons  (fall  versus  spring),  and  from  different  geographic  lo¬ 
cations  (White  Sands  Missile  Range,  Wallops  Island,  and  Fort 
Churchill).  As  can  be  seen  in  Figure  1,  both  the  slopes  and 
amplitudes  of  the  spectra  are  strikingly  similar.  VanZandi 
[1982]  and  Cira  (1984)  were  the  first  to  make  the  suggestion 
that  the  internal  buoyancy  waves  of  the  atmosphere  may  have 
universal  PSDs  and  to  give  the  empirical  basis  behind  it  The 
purpose  of  the  present  report  is  to  present  a  physical  expla¬ 
nation  of  both  the  slope  (n  —  —  3  approximately)  and  univer¬ 
sal  amplitude  of  the  vertical  wave  number  PSDs.  The  theory 
has  an  important,  testable  prediction,  namely,  that  there 
should  be  small  departures  from  a  constant  amplitude  of  the 
PSD,  the  largest  being  that  the  amplitudes  of  the  tropospheric 
PSDs  should  be  smaller  by  a  factor  of  3  or  4  than  those  from 
the  stratosphere.  The  latter  prediction  is  included  in  the  spec¬ 
tral  model  of  Garrett  and  Monk  [1972]:  however,  the  present 
theory  gives  a  physical  mechanism  to  explain  it  as  regards 
atmospheric  spectra. 

2.  Saturation  Model  Due  to 
Kelvin-Helmholtz  Instability 

The  remarkably  constant  PSD  amplitudes  found  in  the 
work  by  Dewan  rt  of.  [1984]  strongly  suggest  that  the  expla¬ 
nation  involves  ‘saturation.'*  In  a  work  by  Phillips  [1977,  sec¬ 
tion  3.4]  there  is  a  description  of  a  saturation  theory  for  the 
wave  PSD.  It  involves  internal  waves  saturating  by  means  of 
shear  or  Kelvin-Helmholtz  instability  induced  turbulence.  The 
first  model  to  be  given  here  is  closely  related  to  Phillips’s 
theory;  however,  it  differs  from  it  in  several  ways.  For  exam¬ 
ple,  we  assume  a  constant  buoyancy  frequency  over  the  alti¬ 
tude  of  interest,  whereat  Phillips  assumes  it  to  be  appreciable 
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only  within  a  narrow  thermodine  and  negligible  outside  of  it. 
In  addition,  we  consider  a  continuum  of  modes,  whereas  he 
considers  only  the  lowest  mode.  Finally,  we  consider  the  one- 
dimensional  vertical  wave  number  spectrum  of  the  horizontal 
velocities,  whereas  he  obtained  the  two-dimensional  horizon¬ 
tal  wave  number  spectrum  for  the  vertical  displacement  of  a 
surface  in  the  thermodine.  Despite  these  differences,  the  key 
dement  (namdy,  saturation  by  means  of  shear-induced  turbu¬ 
lence)  is  exactly  the  same.  In  addition,  our  mathematical  ap¬ 
proach  will  be  identical  to  that  of  Phillips.  Our  second,  but 
related,  model  involves  saturation  through  convective  insta¬ 
bility  induced  turbulence.  Both  models  give  the  same  spectral 
shape  and  dependence  on  buoyancy  frequency.  Convective 
instability-induced  saturation  and  saturation  in  general  have 
been  extensively  reviewed  by  Friits  [1984].  The  first  treatment 
of  saturation  by  convective  instability  was  by  Hodges  [1967]. 
In  section  5  a  more  exact  but  less  physically  obvious  approach 
will  be  given. 

In  the  following  we  shall  very  briefly  review  some  basic 
information  that  we  need  for  this  problem.  The  reader  will 
find  extensive  discussions  by  Rayleigh  [1883],  Gi/f  [1982],  Phil¬ 
lips  [1977,  section  3.2],  Turner  [1973],  and  Ugh  thill  [1978]. 
We  commence  by  listing  the  assumptions  that  we  shall  need. 
These  too  are  discussed  at  length  in  the  cited  references. 
(I)  The  earth's  rotation  has  negligible  effects  [see  GUI,  1982,  p. 
207].  (2)  The  fluid  is  statically  stable.  (3)  The  wave  frequency, 
u,  is  less  than  the  buoyancy  frequency,  N.  Finally,  (4)  we  shall 
assume  the  validity  of  the  Boussinesq  approximation. 

Regarding  assumption  (4),  it  implies  the  following  con¬ 
ditions.  (1)  Effects  of  variations  in  density  are  negligible  with 
regard  to  acceleration  or  inertial  effects.  (2)  The  velocity  fluc¬ 
tuations  are  small  enough  to  make  nonlinear  effects  negligible. 
(3)  Density  fluctuations  consist  of  small  deviations  from  equi¬ 
librium  values.  (4)  The  vertical  scale  of  motion  or  vertical 
wavelength  is  much  less  than  a  scale  height,  H.  We  shall  also 
ignore  viscous  effects  and  assume  that  the  wavelengths  are 
small  enough  to  allow  one  to  ignore  acoustic  effects.  This  Isst 
condition  is  assured  by  the  fact  that  our  wavenumbers,  k,  ere 
significantly  larger  than  N/c,  where  c  is  the  speed  of  sound 
and  N  is  the  buoyancy  frequency.  Also,  we  assume  that  vari¬ 
ations  of  density,  p,  with  respect  to  altitude,  z,  are  much 
smaller  than  variations  of  vertical  velocity  with  respect  to 
height  (e  more  explicit  form  of  condition  (4)  above).  Finally, 
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Fig.  1.  PSD»  from  tropospheric  sad  slralotptwric  nmunanu 
of  vertical  profiles  of  boruoaul  winds:  Dashed  curves  (tropospheric) 
are  from  Donirlx  [I9C2]  and  tMtkk  md  Singleton  [1969]  (sac  Figure 
2.  curve  marked  November);  solid  curves  (stratospheric)  are  from 
Arams  *l  al.  [19S4],  With  the  exception  of  Endlicb's  PSD.  these  PSDs 
are  the  fitted  lines.  The  stratospheric  curves  are  based  on  Table  I  but 
have  been  extrapolated  to  i  *•  100  m. 


we  leave  out  of  account  mean  winds  and  consider  only  Use 
perturbation  velocities. 

Newton's  second  law,  in  terms  of  fluctuations,  is  then  given 

by 


du 

_ar 

(i) 

Poai~ 

dx 

dv 
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p°7,' 

~dy 

dw 
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0) 

where  is  the  equilibrium  density,  p'  is  the  departure  from 
Pa.  0  is  the  acceleration  of  gravity,  r  is  the  departure  of 
pressure  from  equilibrium,  ( is  the  time,  and  «,  e,  and  w  are 
components  of  velocity  fluctuations  in  the  a,  y,  and  i  direc¬ 
tions,  respectively.  The  equations  of  continuity  for  the  incom¬ 
pressible  fluid  are 
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We  aeal  employ  the  standard  procedure  to  etiminale  the 
premise  and  density  terms.  Details  are  to  be  found  in,  for 
eeeeeple,  the  work  by  Glfl  [1912,  p.  139].  One  thus  finds  the 
well  known  equation: 

jp<V*wH.W*V »-0  (6) 


where 


«<Po 

it 


(7) 


is  the  buoyancy  frequency  squared  [Rayleigh,  1883].  With 
compressible  fluids,  a  term  ( -  g/c*)  is  added  to  the  right  side 
of  (7).  The  approximation  is  used  that 


i  a  /  a  \  a*w 
Z7t{p'7tw)*7? 


(8) 


(this  corresponds  to  (kj » |1/W[)  and  where 

0*  a*  a* 
v  "  a** +  ay* +  a*» 


il  £. 

ax* +  ay* 


(9) 


Equation  (6)  is  very  closely  related  to  the  Taylor-Goldstcin 
equation. 

Under  our  assumptions  (N*  is  constant  with  respect  to  «)  it 
is  known  that  (6)  possesses  wavelike  solutions  of  the  form 
[Phillip*,  1977] 


w  —  Wa  exp  {qkgX  +  k,y  +  fc,z  —  on)} 
provided  that  the  following  condition  is  met: 

o£  k ,*  +  k ,* 

N*  “  k.*  +  k,*  +  fc.* 


(10) 

(ID 


where  the  right  side  it  the  square  of  the  cosine  of  the  angle 
between  k  and  the  horizontal. 

As  for  the  horizontal  velocities,  call  them  d*.  they  also  have 
a  wave  solution  of  form 


r,>lg  exp  iHk^x  +  kry  +  kIz  —  wt)} 
provided  (from  (4)) 


k-v-0 


(12) 

(13) 


where  v  is  the  three-dimensional  fluctuation  velocity.  In  other 
words,  particle  motion  is  limited  to  directions  which  arc  per¬ 
pendicular  to  the  wavenumber  vector.  Assuming  conditions 
(13)  and  (II)  and  hence  the  validity  of  (12)  and  (10),  we  can 
write  down  tbe  following  relations  which  are  crucial  for  the 
remainder  of  our  argument: 


— 

it 


Ikjtt 


(14) 


We  now  turn  to  wave  saturation  and  commence  with  the 
one  due  to  shear  instability.  At  the  outset  it  should  be  men¬ 
tioned  that  we  treat  the  shear  instability  here  mainly  for  rea¬ 
sons  of  making  our  explanation  dear.  It  is  easier  to  consider 
the  convective  insubi'ity  (iA,  the  ooe  which  is  clearly  the 
more  probable  one,  since  it  is  4  times  lower  in  threshold)  as  a 
second  step  in  the  discussion.  As  has  been  pointed  out  orig¬ 
inally  by  Miles  [1961,  1963]  and  Howard  [1961],  tbe  Rich¬ 
ardson  number,  Hi  (defined  at  W*/S*.  where  S  is  tbe  vertical 
shear  of  tbe  horizontal  velocity),  obeys  the  following  relation 
as  tbe  necessary  condition  for  dynamic  instability  (leading  to 
turbulence): 


*.<*  (IS) 

We  shall  follow  FMfllps  [1977]  and  take  this  as  being  also  the 
sufficient  condition  for  instability.  This  is  valid  for  practical 
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purpoM*  evaa  though  il  is  sot  rigorously  true  mathematically. 
Thus  the  condition  on  (hear  at  the  border  line  of  turbulence  i» 


Note  that  the  subscript  "mu’*  iinpiie*  that  (16)  hold*  for  the 
maximum  shear*  along  a  wave  train.  Equation  (16)  can  be 
simplified,  without  lot*  of  generality,  by  choosing  a  coordinate 
system  so  that  one  horizontal  axis,  say  the  y  axis,  is  in  the 
direction  of  the  horizontal  velocity.  It  can  be  written  then  as 

In  order  to  closely  examine  the  physical  significance  of  (17) 
we  consider  a  numerical  example.  Let 

n  -  n»  sin  M  («) 

At  critical  shear,  (IS)  gives 


[&,Uo  cos  k 


-I 


2 N 


(19) 


where  the  argument  of  the  cosine  must  be  O’  or  ISO*.  For 
specificity,  let  N  «•  2  x  10**  s'*  and  2,  -  500  m.  From  (19) 
the  critical  velocity  amplitude  is  defined  as 

(-Jl»*^(2N)- 3.18  m/s  (20) 

If  u,  had  exactly  the  value  3.18  m/s,  the  wave  would  be  at.  but 
not  over,  the  threshold  of  instability.  In  order  that  the  insta¬ 
bility  grow  in  a  finite  time,  the  value  of  u«  must  exceed  3.18 
m/a  Consider  the  case  of  «,«  3.19  m/s.  It  is  easy  to  show 
from  (19)  that  the  shear  will  exceed  the  critical  value  over  two 
small  regions,  of  depth  10  m,  in  a  wavelength  of  900  m.  In 
other  words,  each  turbulent  region  in  a  wavelength  will 
occupy  2%  of  the  wavelength,  and  there  will  be  only  two  such 
regions  in  one  wavelength.  From  this  illustration  it  is  dear 
that  (1)  the  thick  nest  of  a  turbulent  layer  can  be  very  much 
smaller  than  the  wavelength  of  the  wave  that  caused  it  (this  is 
true  despite  the  fact  that  subsequent  spreading  or  other  thick¬ 
ening  mechanisms  are  possible)  and  (2)  there  is  no  fixed  rela¬ 
tion  between  2  and  turbulent  layer  thickness.  Finally,  it  is 
dear  that  waves  which  have  velocity  amplitudes  above  thresh¬ 
old  will  lose  energy  to  turbulent  dissipation  until  they  reach 
the  critical  limit.  When  the  latter  occurs,  the  dissipation  mech¬ 
anism  will  no  longer  operate.  This  is  a  classic  case  of  “satu¬ 
ration.’* 

Our  goal  it  to  relate  the  one-dimensional  PSD,  6  (k,)  to  the 
saturation  condition.  First,  we  wish  to  start  (from  (20))  with 


<«eLs 


(21) 


Aa  mentioned  above,  we  duplicate  the  procedure  given  by 
Phillips.  The  first  step  is  to  review  the  fact  that  (as  is  well 
known)  the  total  variance  of  the  fluctuations  is  equal  to  the 
PSD  intagratnd  over  all  wave  numbers  [Phillips,  1977,  section 
3.4,  p  220). 

f  flfcjdk.-  I<uHkJ>kr-4«<N*>  (22) 
JUt,  • as. 

In  the  special  case  o fa  aiae  wave  (e*.  (180  which  has  a  single 
wave  anther  it  is  wen  known  that 


**,*•<«*>  (23) 


Of  course,  if  one  considers  a  band  of  Fourier  components 
around  a  given  fc,  the  variance  <u*>  would  not  be  exactly 
equal  to  of  the  main  component  However,  the  variance 
of  such  a  hand  would  be  proportional  to  u0*, 

(24) 


For  such  a  band  the  integral  of  the  PSD  reduces  to 


u 


WJM. 


(25) 


Next,  we  again  follow  Phillips  and  assume  that  the  band 
width,  6kv  is  proportional  to  the  central  wave  number,  in., 


Alt,  oc  it,  (26) 

and  we  turn  our  attention  to  the  case  where  there  are  many 
waves  containing  such  bandwidth*.  As  Phillips  put  it,  one  can 
visualize  the  field  of  gravity  waves  as  “a  random  succession  or 
wave  groups  at  different  wavenumbers.”  In  this  case,  the  con¬ 
tribution  to  the  mean  square  velocity  fluctuations  from  each 
wave  number  region  is 

0 (27) 


where 

<u*>  ■  £  <u»(kj>w^-  -  f  flfcj  dk,  (28) 
uu.  JUu. 

Neat,  as  in  the  work  by  PkiUips  [1977],  we  make  the  assump¬ 
tion  that  all  the  wave  number  components  are  limited  individ¬ 
ually  by  the  saturation  condition,  (17).  We  then  have,  with  the 
help  of  (24)  and  (27X 

(29) 


and,  from  (21). 


4N* 

*.* 


(30) 


This  is  the  result  we  were  seeking. 

Note  that  in  (30)  one  bat  a  proportionality,  not  an  equality 
(as  waa  also  true  in  Phillips’s  treatment).  Note  also  the  basic 
assumptions:  (1)  each  wave  number  band  saturates  by  the 
shear  instability.  (2)  the  wave  hands  dk,  scale  as  k„  and  (3)  the 
wave  bands  saturate  independently.  These  assumptions  of 
Phillips  have  met  with  success  for  the  case  he  studied.  How¬ 
ever,  them  assumptions  arc  not  needed  in  the  more  exact 
treatment  shown  in  section  3  below. 

There  is  a  for  shorter  and  more  direct  method  to  obtain 
(30).  It  it  by  the  “principle  of  similitude"  or  “dimensional 
analysis.”  This  approach’s  validity  rests  completely  on  the  ini¬ 
tial  assumption  made.  From  (15)  we  can  assume  that  if  satu¬ 
ration  due  to  the  shear  instability  determines  p,  then  p  must 
depend  upon  the  critical  shear  S,1  —  4 N1  (let  N  be  constant). 
The  only  other  variable  that  p  would  depend  upon  in  this  case 
is  kr  Using  the  fact  that  the  dimensions  of  p  are  [L*T~  *]  and 
those  of  S,  and  k,  art  [T*1]  and  [f*],  respectively,  one 
immediately  obtains 


P* 


S* 

k.* 


4 N* 

V 


(31) 


where  a  universal  constant  of  proportionality  of  order  unity  is 
left  to  be  determined  tty  means  of  experiment.  This  does  raise 
the  question  of  how  many  of  the  assumptions  made  above  am 
actually  necessary  fo  arrive  at  (31).  On  the  other  hand,  it 
appears  that  our  urn  of  Ak,  oc  k,  it  a  crucial  one,  and  the 
dhneneiowal  argument  leaves  this  rather  well  hidden.  It  te  for 
this  reason  that  we  have  given  both  derivations  in  this  paper. 


2745 


_ vf  and  Qooo:  ATtsosrttsatc  Scalar  Winds  -^y 


Fig.  2.  After  Endfick  and  Singleton  [1969,  Figure  10].  Smoothed 
•poem  for  mmcB  of  wind  profile*  on  April  I,  1966,  July  5,  1966. 
•ad  November  10,  1965.  The  “November”  curve  it  plotted  to  2  —  5 
km  ia  Figure  I  for  comparison, 


Note  that  I  umlty  [1964]  arrived  at  a  similar  dependence  to 
that  given  in  (31),  for  d.  but  he  did  to  in  an  entirely  different 
context,  namely,  that  of  buoyancy  subrange  turbulence.  All 
dimensional  arguments  are  subject  to  possible  omission  of  un¬ 
derlying  physics,  as  this  case  illustrates. 


1  Saturation  Model  Due  to 
Convective  Instability 

Hodges  [1967]  was  the  first  person  to  suggest  that  the  satu¬ 
ration  due  to  convective  instability  could  limit  the  amplitudes 
of  upward  propagating  gravity  waves.  Frit  it  [1984]  has  given 
an  excellent  review  of  this  work  and  related  subjects.  We  now 
consider  the  effect  of  this  particular  mechanism  upon  jikj,  the 
vertical  wave  number  PSD. 

The  convective  instability  condition  [see  Frills,  1984]  is 
given  by 


di  iff  „ 
dr  dx 


(32) 


where  6  and  ff  are  the  mean  and  perturbation  values  of  the 
potential  temperature.  Using  a  truncated  Taylor  expansion,  it 
is  easy  to  show 


9  - 


(33) 


and  using  (32)  with  the  equality  sign  to  denote  the  onset  of 
instability,  we  obtain  the  critical  velocity  as 

<wL  -  “  (37) 

To  simplify  calculations,  let  us  consider  this  in  the  two  dimen¬ 
sions  of  x  and  z,  Le,  let  u  stand  for  horizontal  velocity.  Then, 
from  (13)  and  (37), 

-MwL-MJt.  (38) 

Using  (38)  in  (37)  gives 

(«L-f  (39) 

Le.,  u  is  equal  to  the  horizontal  phase  velocity,  and  using  (11) 
with  it,  —  0, 

'“'•-iTTT?3'©  m 

where  we  have  followed  the  argument  of  Frills  [1984]  that 
k, »  k. 

At  this  point  one  needs  merely  to  repeat  the  arguments 
following  (21)  in  order  to  arrive  at 

iHkJoc^s  (41) 

which  shows  that  the  saturation  amplitude  is  4  times  smaller 

than  would  be  given  by  the  shear  instability  mechanism.  We 
shall  assume  that  the  convective  instability  controls  saturation 
in  the  more  exact  treatment  of  section  S  below.  (See  Cossard 
and  Hooke  [1973]  for  further  discussion  of  this  instability.) 
The  approach  for  the  shear  inilabiity  would  be  altered  by 
replacing  4AI*  by  Na.  In  addition,  the  same  dimensional 
analysis  derivation  would  be  valid,  as  was  given  above.  There 
seems  to  be  little  question  that  the  convective  mechanism  is 
the  controlling  factor  in  saturation  when  compared  to  the 
shear  mechanism. 

4.  Comtabison  With  Extexi  mental  Results 
4.1.  PSD  Log-Log  Slopes 

Endttch  and  Singleton  [1969]  published  PSDs  from  vertical 
profiles  of  horizontal  winds  that  were  measured  by  means  of 
radar-tracked  jimspbercs.  These  measurements  were  in  the  200 
m  to  16  km  altitude  region,  and  the  PSDs  covered  a  wave¬ 
length  range  of  100  m  to  20  km  (see  Figures  1  and  2).  Over  the 
wavelength  range  of  100  m  to  5  km  they  observed  slopes  of 
—  }.  When  the  range  is  reduced  to  100  m  to  7  km,  the  slope  is 
ctoac  to  —3,  in  agreement  with  our  theory.  At  very  large 
wavelengths  (2  >  5  km)  they  saw  a  pronounced  change  (sec 


where  k  is  the  vertical  displacement  of  the  fluid  parcel  causing 
the  fluctuation  9.  but,  with  the  help  of  (10), 

<«! 

hence,  from  (33) 


'-©(=) 

Taking  the  derivative  of  (33)  with 


d»  dt/ik.\ 
dtmdt[ko)W 


(35) 

to  z  with  dt/dx 

(36) 


TABLE  1.  Fitted  Log  PSD  Amplitudes  and  Slopes  of  PSDs  From 
Dtma a  H  aL  [1984]. 


Label 

Wiad  Profile 

Log  PSD  Amplitudes 
(From  Fitted  Line) 
at  2  ■>  I  km, 

(m/s)Mc/at) 

Slope 

A 

May  22,  197t 

340 

—  349  ±  0.16 

B 

May  24  1978 

3.28 

-241  ±  021 

C 

April  26,  1977 

3.41 

—2.99  ±  009 

D 

May  2. 1971 

327 

—  2.91  ±  OI9 

E 

Sept.  12, 1978 

3.55 

—  326  t  007 

Average 

3,42  ±  —0.15 

-3.0  ±  0.33 

The  raagt  is  from  1>  I  km  to  2  •>  200  m. 
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TABLE  2.  N*  Scaling  of  Spectral  Amplitudes 


Altitude 

Bangs, 

km 

Region 

A*  x  10-4. 
(rad/s)1 

A. 

r*d/» 

Ratio  of 
A.VA1 
Upper 
Suatospherc 

Inverse 

Ratio 

0-10 

Troposphere  (),  6,  10  km) 

132 

onus 

0286 

330 

12-20 

Lower  stratosphere 

441 

00210 

0955 

IBS 

2I-S0 

Upper  stratosphere 

462 

00215 

1 

1.00 

50-72 

Lower  Tnnosphfft 

274 

00166 

0593 

1.69 

72-86 

tipper  onoiptoi 

3j63 

0786 

1372 

86-91 

Thermosphere 

486 

lid 

0952 

Baled  oa  11$  Standard  Aunoapbere  (1976). 

Figure  2).  It  is  accessary  to  choose  wavelengths  less  than  1  km 
to  eliminate  effects  bom  the  long  wavelengths  in  order  to 
observe  the  saturation  effects. 

Devon  et  al.  [1984]  reported  PSDs  of  such  wind  profiles 
measured  by  the  smoke  bail  method  in  the  altitude  range  of 
13-37  km.  They  reported  a  slope  of  —3.07  ±  0.3  in  the  wave¬ 
length  range  of  300  m  to  1  km,  in  agreement  with  theory  (see 
Table  1).  They  also  reported  a  slope,  if  a  line  were  fitted  bom 
1  km  to  40  m.  of  —2.7  ±  0l2;  however,  in  the  upper  end  of 
their  spectra,  which  had  a  Nyquist  wavelength  of  30  m, 
showed  much  curvature.  By  fitting  up  to  40  m,  the  half  Ny- 
quiit  value,  it  was  presumed  that  aliasing  effects  would  be 
negligible;  however,  other  artifacts  such  as  measurement  error 
noise  could  extend  to  larger  wavelengths,  (n  view  of  this,  one 
may  possibly  regard  the  slope  over  1  km  to  200  m  as  being  the 
more  valid  slope.  On  the  other  hand,  sac  section  5,  which 
shows  that  the  saturation  hypothesis  is  compatible  with  both 
of  these  slopes. 

Rosenberg  ft  aL  [1974]  also  reported  PSDs  baaed  on  NASA 
smoke  trail  measurements  (vertical)  in  the  altitude  range  of 
5-18  km  and  wavelength  range  100  m  to  3.6  km.  They  cited 
an  average  slops  of  — Z7S  ±  Oil.  If  they  had  restricted  their 
wavelength  range  to  the  interval  I  km-100  m,  then  one  would 
ffHpfct  a  looewhat  steeper  slope.  This  decrease  in 

slope  at  smaller  wave  number  can  be  seen  in  Figure  2,  which 
is  based  on  Endlich's  data. 

Finally.  Daniels  [1982]  reported  on  a  total  of  1200  vertical 
profile  measurements  made  with  rawinsondes  in  the  4-  to 
16- km  altitude  range  and  wavelengths  between  4  km  and  80 
m.  His  spectral  slopes  were  in  the  neighborhood  of  —2-5  (see 
Figure  1).  Since  only  the  resulting  fit  and  not  the  actual  PSD 
is  shown  by  Daniels,  we  believe  the  determination  of  the  slope 
was  dominated  by  the  long-wavelength  region.  If  the  wave¬ 
length  range  were  hmited  to  I  km  to  80  m.  then  one  would 
expect  the  slops  to  approach  —3,  as  was  the  case  with  the 
Endlich  data.  Either  slope,  as  is  shown  in  section  S,  is  compat¬ 
ible  with  saturation. 

In  conclusion,  observations  of  spectral  slope  over  the  100  m 
to  1  km  range  appear  to  be  consistent  with  the  theoretical 
value  of  —3  in  the  sense  that  when  all  the  observations  are 
carefully  considered,  this  slope  is  either  directly  supported  or, 
at  least,  not  ruled  out.  Contamination  bom  long  wavelength 
power  it  a  problem,  and  cartful  PSD  analysis  is  required  to 
arrive  at  the  theoretical  slope. 

42.  PSD  Amplitude! 

Aswan  [1964]  rapened  the  spectral  slopes  fitted  to  the  PSD 
la  the  range  of  1  km  in  200  at  The  average  of  the  log  ampli¬ 
tude  of  the  feted  fine  (to  the  PSD)  at  2  «  1  km  is  3.42  ±  0.15 
(sse  Table  1).  This  small  variability  is  good  evidence  for  satu¬ 
ration  indeed.  Equation  (30)  however,  implies  that  +  (k J  scales 
as  N\  which  is  to  say  that  the  amplitude  is  only  approxi¬ 


mately  constant  in  that  it  should  vary  between  altitude  regions 
where  N *  differs  significantly.  Table  2  lists  values  of  A*  for 
different  altitude  ranges.  As  can  be  seen,  a  significant  differ¬ 
ence  in  A*,  and  hence  p  amplitude,  should  exist  between 
stratospheric  and  tropospheric  values.  Theory  predicts  a 
(actor  of  3.5  difference  in  that  case,  the  troposphere  having  the 

fyriilW 

As  can  be  seen  bom  Figure  1,  the  value  of  Daniels's  [1982] 
PSD  at  wavelength  equal  to  I  km  (partly  in  the  troposphere) 
is  indeed  significantly  smaller  in  relative  value  than  the  strato¬ 
spheric  values  of  flows  et  al.  [1984].  Endlich  and  Singleton's 
[1969]  value  (alls  between  them.  A  cross  it  located  at  theoreti¬ 
cal  value  using  Anor*  ”  1-32  x  10" 4  red*  »"*.  The  theoreti¬ 
cal  tropospheric  average  (which  is  calculated  from  the  strato¬ 
sphere  average  by  multiplying  the  latter  by  hfnar*/h/mAT*)  “ 
in  perfect  agreement  with  experiment.  This  observation  can  be 
considered  as  strong  evidence  that  the  PSD  is  determined  by 
saturation. 

Finally,  it  should  be  mentioned  that  the  strongest  evidence 
for  our  general  saturation  hypothesis  lies  in  the  following  ob¬ 
servation.  The  amplitude  of  upward  propagating  buoyancy 
waves  if  unsaturafed  increases  at  the  inverse  of  the  square  root 
of  the  atmospheric  density  (p),rt  (see  Costard  and  Hooke 
[1975,  p.  76].  This  fact  implies  that  the  PSD  amplitude,  for 
any  given  wave  number,  would  increase  with  altitude  as  the 
inverse  of  p.  As  can  be  seen  bom  Tabic  3.  this  would  cause 
PSD  amplitudes  to  increase  rapidly  as  a  function  of  altitude. 
Sines  the  experimental  results  of  Figure  t  rule  out  such  a 
dependence,  tbs  saturation  hypothesis  is  strongly  supported. 
(Note,  however,  that  the  smoke  trail  method  used  for  velocity 
PSD’s  was  Hmited  to  fair  weather  conditions.) 

5.  The  Bxoad-Band  Saturation  Model 

As  has  been  previously  mentioned,  the  “narrow  band”  treat¬ 
ment  given  above  (wherein  each  band  is  assumed  to  saturate 
independently)  is  certainly  only  an  approximation.  The  pur¬ 
pose  of  this  section  is  to  take  into  account  the  effects  of  all  the 
other  bands  upon  the  saturation  of  a  given  hand. 

In  order  to  ascertain  the  mean  square  shear  due  to  the 
fluctuations  one  can,  by  definition,  integrate  the  shear  spcc- 


TABLE  3.  Density,  *  FsBoff  With  Altitude 


Altitude, 

km 

ftegioa 

P* . 
(kg/m*) 

WAreos)' 1 

0 

Troposphere 

133 

t 

10 

Lower  stratosphere 

0414 

197 

30 

Mid  stratosphere 

1J4  x  10** 

66.9 

50 

Stratosphere 

1433  x  10*» 

1.19  x  I0» 

90 

nenMNpbcre 

342  x  I0*» 

33  x  10* 

•US  Standard  Atmosphere  (1976). 


Dew  \m>  Good:  Aneoanoauc  Scalar  Wind* 


2747 


TABLE  A  Sensitivity  Test  for  Model 


Variable 

k,,  i/m 

c/m 

atm'/t'Hc/m) 

m 

Comments 

*. 

... 

15  x  10*» 

3  x  I04 

... 

a  iacmsMS  a  (actor 

2  x  10*' 

IS  x  10-* 

5  x  10* 

146 

of  2  far  a  focior  of 

3  x  10** 

IS  x  10*' 

5  x  10* 

34)7 

3  increase  in  k. 

1  x  10*4 

IS  x  tO'* 

5  x  10» 

5.97 

(quite  sensitive) 

*— 

5  x  to-' 

... 

5  x  10T* 

... 

a  increases  a  (actor 

3  x  10*' 

I0‘» 

5  x  10T4 

168 

of  Ml  for  3 

3  x  10"' 

3  x  10-* 

5  x  107* 

196 

orders  of  atagoitude 

3  x  10*' 

5  x  10-* 

5  x  10f* 

3.10 

increase  in  k_„ 

5  x  tO*' 

to-* 

5  x  10^* 

112 

(not  sensitive 

3  x  10'* 

uo 

5  x  ur* 

3.16 

for  this  range  of  a) 

B 

3  x  10*' 

15  x  10-* 

... 

... 

a  increases  a  (actor 

3  x  tO’' 

13  x  10*' 

5  x  10* 

145 

of  1.3  for  a  factor 

3  x  I0-» 

13  x  10*» 

1  x  10* 

161 

of  16  increase  in  B 

5  x  10-' 

13  x  10"* 

2  x  10* 

179 

(mild  sensitivity) 

5  x  tO-' 

15  x  10*' 

3  x  I04 

190 

5  x  tO'* 

15  x  10-' 

1  x  10* 

327 

N‘  -  3.35  x  Hr*  Hz1  for  All  cam.  Time  results  can  ba  undemood  Bom  /,*-***  4k  -  N‘.  A 

steepened  dope  (larger  a)  decteaiei  the  integral,  while  the  increaae  of  k _ and  B  increases  the 

integrals 


tram,  fc*p,  over  all  wave  numbers.  The  borderline  of  the  con¬ 
vective  instability  occurs  (under  the  assumptions  stated  pre¬ 
viously)  when  JU  ■»  1.  Hence  one  arrives  at  (within  a  constant 
very  dose  to  unity) 

f*~  fc'*  dk-N1  (42) 

The  upper  limit,  k _ is  the  wave  number  where  the  waves  has 

been  replaced,  for  the  most  part,  by  turbulence 

There  is  evidence  that  for  wave  numbers  less  than  a  certain 
value,  k.,  the  spectrum  bends  over  from  the  constant  slope 
regime  to  a  regime  which  is  flat  and  horizontal  for  the  oceanic 
case  [Carreti  and  Munk,  1972.  1973],  For  the  atmospheric 
case,  T.  E.  VanZandt  (private  communication,  1983)  has  sug¬ 
gested  the  analytical  expression  of  the  form 

k-i'+B'  <«> 

In  the  following  we  shall  understand  that  the  region  where  the 
PSD  goes  as  k~m  is  primarily  governed  by  saturation,  whereas 
in  the  flat,  horizontal  region  there  are  only  indirect  contri¬ 
butions  to  saturation.  The  shape  of  the  latter  region  is  pre¬ 
sumably  due  to  properties  of  the  wave  sources  as  well  as  to 
mode  interaction,  hut  at  present  this  remains  an  open 
question. 

It  is  sufficient  for  our  purposes  to  represent  this  situation  by 
means  of  the  simple  expression 

*-£  (0  Zk<kJ 

(44) 

(k.sksfk^j 

where  £  m  Ak,“*.  Eliminating  A  and  inserting  (44)  into  (42) 
and  integrating,  one  obtains 

+  {*«,»-- VJ  (43) 

where  a  *  3  but  where  a  can  be  very  dose  to  3  so  for  as 
numerical  calculations  are  concerned. 

Equation  (43)  can  be  easily  solved  numerically  for  a  for 
given  values  of  N\  £,  k,.  and  kmm.  In  particular,  ws  determine 
a  from  parameters  based  on  Figure  2.  There  we  iafor  that 
k,  •  3  x  10**  c/m  (cydes/meter)  and  i»5*  10* 
(mVc’Vte/m).  approsunatdy.  One  can  obtain  the  average 
value  of  N  for  the  troposphere  from  Table  2;  however,  this 
value  must  be  divided  by  2a  to  make  the  units  compatible 


with  those  of  it  (Le,  c/m  as  opposed  to  rad/m);  hence  we  take 
N*  -  3.33  x  10-*  Hz*.  For  we  shall  use  the  value  (1/40 
m)  on  the  basis  of  Demon  et  at.  [1984],  While  it  is  true  that  the 
latter  measurements  were  carried  out  in  the  stratosphere,  we 
dial!  assume  that  the  troposphere  does  not  differ  significantly 
in  this  regard.  As  we  shall  see  in  Table  4,  the  results  for  x  %  3 
arc  very  insensitive  to  the  value  of  k _ hence  this  approxi¬ 

mation  is  not  problematic.  Solving  for  a,  we  found  n  —  3.07. 
From  this  result  we  find  good  evidence  that  in  regard  to  the 
value  of  n  the  narrow-band  approximation  is  a  valid  one. 

Table  4,  as  has  already  been  mentioned,  explores  the  sensi¬ 
tivity  of  n  in  (43)  to  changes  in  the  input  parameters.  While 

there  is  little  sensitivity  to  k _ there  is  moderate  sensitivity  to 

£  (or,  equivalently,  to  N2)  and  a  large  sensitivity  to  k.. 

In  regard  to  the  high  sensitivity  to  k„  values  we  must  re¬ 
spond  to  a  remark  made  to  us  by  D.  Fritts  in  private  conver¬ 
sation  (1983).  He  suggested  that  at  low  values  of  k  in  Figure  2 
there  could  be  significant  contributions  from  nongravitational 
wave  sources  such  as  the  jet  stream.  For  this  reason,  be  sug¬ 
gested  that  one  estimate  1,  and  £  from  an  alternate  figure  in 
Endiich  and  Singietotft  [1969,  Figure  3]  paper,  in  which  the 
PSD  is  calculated  on  the  basis  that  mean  winds  (ie,  the 
average  over  a  number  of  profiles  taken  on  different  days)  are 
removed  by  subtraction.  On  this  basis  we  obtain  £ »  10* 
(m*/sa)/Hz  and  k.  -  3  x  10‘4  (c/m).  If  we  retain  k_„  -  (1/40 
m)  and  N*  -  3.35  x  10‘*  Hz’,  we  arrive  at  n  -  2.70.  This 
value,  while  lower  than  3.0,  is  not  very  much  lower,  and  it  is 
still  reasonably  consistent  with  slopes  observed  by  Endiich 
and  others  [e.g*  Daniels,  1982]  in  the  troposphere. 

As  we  see,  the  present  and  more  accurate  model  leads  to 
values  of  a  that  are  not  exactly  equal  to  3  but  which  are 
nevertheless  consistent  with  our  saturation  hypothesis.  With 
this  in  mind,  the  values  of  a  -  2.5  and  2.73  ±  .1  reported  by 
Daniels  [1982]  and  Rosenberg  et  of.  [1974],  respectively,  etc. 
may  possibly  not  be  due  merely  to  artifact  as  was  suggested 
previously.  They  may  possibly  have  actually  been  closer  to  the 
slopes  “in  nature.”  In  any  case,  it  is  important  to  note  that  the 
broadband  modd  permits  values  of  a  which  are  different  from 
3  but  that,  to  apply  this  model,  one  must  have  information 
about  k.,  k__  £,  and  N, 

In  passing  it  should  be  mentioned  that  the  spectrum  found 
by  Garrett  and  Monk  [1972,  1973]  for  the  upper  ocean  (with, 
in  their  case,  a  2.3)  may  be  entirely  consistent  with  our 
explanation.  Their  PSD  amplitudes  certainly  scale  as  N*.  We 
intend  to  treat  this  in  detail  elsewhere  in  the  future. 
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Out  saturation  hypothesis  suggests  other  testable  predic- 
titMit.  Some  of  these  have  been  pointed  out  to  us  by  Smith, 
Fritta,  and  VanZandt,  who  plan  to  publish  them  in  the  near 


rWnirtt.  G„  Terrestrial  ant 


(climatic)  eritsris  guidelines  for 
I  M2  Revision.  NASA  Tech. 


future  in  Geophysical  Research  Let  lets. 


6.  Summasy 

A  simple  physical  explanation  for  the  nearly  universal  verti¬ 
cal  wave  number  PSD  for  horizontal  winds  has  been  present¬ 
ed.  It  rests  on  the  plausible  assumption  that  such  PSD  ampli¬ 
tude*  and  shapes  arc  determined  by  saturation  due  to  turbu¬ 
lent  breakdown.  As  was  pointed  out,  this  general  approach 
has  been  used  in  the  physics  of  ocean  surface  waves  and  ther- 
modine  interna]  waves.  The  theory  predicts  a  PSD  slope  in 
the  neighborhood  of  -11  to  —3  and  an  amplitude  depen¬ 
dence  that  scales  as  N*.  These  predictions  are  in  agreement 
with  experiment  The  N1  dependence  as  well  u  the  approxi¬ 
mate  -3  slope  gives  us  two  ways  to  test  the  theory.  In  partic¬ 
ular,  the  predicted  difference  of  3  to  4  between  tropospheric 
and  stratospheric  PSD  amplitudes  seems  to  provide  a  practi¬ 
cal  test 

The  simplicity  of  our  model  is  one  of  its  most  important 
assets.  It  may  also  provide  a  physical  explanation  of  tbs  verti¬ 
cal  wave  number  PSD  of  the  universal  spectral  model  pro¬ 
posed  by  Garrett  and  hitmk  [1972,  1975]  for  mesoscale  ocean¬ 
ic  motion*,  at  least  lor  those  situations  and  wave  numbers 
where  saturation  can  be  proved  to  be  a  dominant  factor. 

In  the  simple  initial  approach,  our  main  assumptions  in¬ 
cluded  (I)  ktotk„  (2)  AA,  qc  *„  and  (3)  individual  Fourier 
wave  components  saturate.  Our  more  exact  approach  did 
away  with  the  need  of  these  assumptions,  and  we  found  that  s 
can  he  different  from  3.  Future  publications  will  address  fur¬ 
ther  predictions  of  our  approach. 

In  this  paper  our  attention  has  been  focused  on  a  range  of 
wavelengths  that  included  40  m  to  about  3  km.  At  the  small 
scale  cod  it  is  assumed  that  turbulence  effects  become  impor* 
tant  (at  A_..l 
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